Quantum Mechanical Fluctuations at the End of Inflation 
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During the inflationary phase of the early universe, quantum fluctuations in the vacuum generate 
particles as they stretch beyond the Hubble length. These fluctuations are thought to result in the 
density fluctuations and gravitational radiation that we can try to observe today. It is possible to 
calculate the quantum-mechanical evolution of these fluctuations during inflation and the subsequent 
expansion of the universe until the present day. The present calculation of this evolution directly 
exposes the particle creation during accelerated expansion and while a fluctuation is larger than the 
Hubble length. Because all fluctuations regardless of their scale today began as the vacuum state in 
the early universe, the current quantum mechanical state of fluctuations is correlated on different 
scales and in different directions. 



I. INTRODUCTION 

When one thinks of macroscopic manifestations of 
quantum mechanics one usually thinks of phenomena 
such as superconductivity and lasers that have crept at 
least somevifhat into our everyday lives. However, if there 
was an inflationary epoch in the early universe, quantum 
mechanical fluctuations and correlations determined the 
large-scale structure of the universe, and galaxies and su- 
perclusters are quantum mechanics writ upon the largest 
scales of our Universe. 

During inflation a slowly rolling scalar field (the infla- 
ton) drives an exponential expansion of spacetime. This 
rapid expansion addresses several cosmological issues. It 
ensures that space is very close to flat and that topo- 
logical defects that form before and early during infla- 
tion will be exceedingly rare today. Finally the quantum 
mechanical evolution of the perturbations in a spatially 
uniform inflaton field (or other fields) provide the seeds 
for structure formation in the recent universe and for the 
fluctuations of the microwave background. 

There are several equivalent ways to calculate the evo- 
lution of fluctuations in scalar fields during the approx- 
imately de Sitter inflationary period. Possibly the most 
straightforward way is to exploit the fact that the equa- 
tion describing the evolution of the flcld operator is the 
same as that describing the classical evolution of a scalar 
fleld fl*]; one can start with a reasonable approximation 
for the wavefunction of the de Sitter vacuum and evolve 
it forward through the de Sitter stage. At the onset of 
radiation domination, one can use the method of Bogol- 
ubov coefficients [2| to ensure continuity of the field and 
its time derivative. 

Studies of how fiuctuations in the early universe lose 
coherence have followed two main paths. The first looks 
at a quantum field coupled to the environment: for ex- 
ample, a thermal bath [3j] or the backreaction of pair 
production of the quantum state of the universe Q . The 
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second path is to examine how the appearance of de- 
coherence develops in de Sitter space without any de- 
cohering interaction. Work to understand how initially 
quantum fluctuations begin to act classically as they are 
stretched outside the Hubble length moved forward dra- 
matically with the work of Guth and Pi [5] . They found 
that long after a mode passes through the Hubble length 
it follows a classical Gaussian probability distribution. 
Several people have addressed the quantum-to-classical 
transition in terms of squeezed states li^, |8| . Albrecht 
et al. [1] argued that although the the concept of squeezed 
states is useful in understanding the evolution of density 
fluctuations, it is entirely equivalent to more traditional 
methods and does not provide any new physical conse- 
quences. 

This paper presents an alternative method that explic- 
itly follows the particle creation. It is perhaps closest in 
spirit to the recent work of Mijic @ that follows the evo- 
lution of particle number as a mode is stretched beyond 
the Hubble length. Rather than treat the evolution of the 
field in the standard manner using a field operator that 
follows the de Sitter vacuum, we can use a operator that 
allows the wavefunction is expanded as a superposition 
of multi-particle states as observed by a small comov- 
ing detector. The entire wavefunction is calculated as a 
function of time, rather than the expectation value of a 
particular operator (typically the field operator). This 
method is computationally intensive, so it cannot be ap- 
plied to modes that spend many Hubble times larger than 
the Hubble length. However, it offers new insights into 
how density fluctuations approach the classical limit. 

In particular one can calculate the particle number us- 
ing both techniques and they agree, but one can also con- 
struct any operator after the wavefunction has been cal- 
culated and determine its expectation value or the prob- 
abilities of the various outcomes of a particular measure- 
ment. [|TT]will derive the Hamiltonian of a scalar field in 
an homogeneous, isotropic universe and define the field 
operators for the de Sitter vacuum ( ijll A[) and the vac- 
uum of a radiation-dominated universe ( jjllSp . The evo- 
lution of the wavefunction is described in i jll CI and its 
initial conditions in i jll Dl 
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The results of the calculation (S|TTT| fall in two cate- 
gories. First, we examine the robust of the calculations 
and compare with the results from the standard treat- 
ments in mil Al We find that this treatment predicts the 
same amount of particle production as in the standard 
treatment ( §11 Al and e.g. [13]) and that the number of 
particles follows a thermal distribution. This highly ex- 
cited Fock state behaves quasi-classically; however, the 
state of the field is still pure and furthermore the phases 
of the various multiparticle states are correlated f ^IIIBp . 
§Vl speculates on the possibility of observing these corre- 
lations and future work. 



II. THE EVOLUTION OF SCALAR-FIELD 
MODES 

We will follow the evolution of perturbations of a uni- 
form real scalar field in a homogeneous and isotropic 
spacetime with flat spatial sections. The metric is 
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where to get the final result we have dropped a term in 
the integrand equal to a total derivative with respect to 
the conformal time. 

The field u(x, r) can be expressed in terms of Fourier 
modes, 



(1) 



where a(T) is the scale factor. The action for the scalar 
field ((/)) is given by 
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where we have taken ?i = c = 1 and \/—g = o-'^{t). We ne- 
glect the back reaction of the perturbations on the metric. 
From the action we would like to find the Hamiltonian 
and determine the quantum mechanical evolution of the 
fluctuations. 

Now let's take = \ (to^ + ^R) 02 + where the 
scalar field may have a non-minimal coupling to gravity 
through the Ricci scalar R [11]. With these we obtain 
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where we have dropped the constant term Vo- Let's make 
the substitution u = acj) to try to absorb the factor of 
a2(r) into the fields. 
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where u-k(''") = "k(''") ensure that the field u is real. 
This yields 



\ I drd^k 



duk 



dr 



(k^ + mls) ]uk] 



(10) 



where k is the comoving momentum of a mode. 

Now we have the action for a related scalar field u in a 
flat space-time with a negative contribution to the square 
of the mass of the mode P, [l^ . 

The value of the mass depends on the evolution of the 
background spacetime (i.e. a(r)). If the dominant energy 
density in the universe is characterised by a pressure that 
is proportional to the energy density {P = wp), The scale 
factor evolves as 
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where H is the Hubble parameter, {da/ dr) / a^ . For w — 
— 1/3, the product of the scale factor and the Hubble 
parameter is constant. The effective mass is 
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The Hamiltonian of the field is given by 



9ul 



dr 
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To look at the quantum mechanics of this Hamiltonian 
one can take one of several routes, depending on the def- 
inition of the field operator Uk. In general it is [l| 



.g(fc,r)ak + g*{k,T)a}_ 
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where and Ok are creation and annihilation operators 
that satisfy the following standard commutator relations 

[4, 4,] = [ak, flk'] = 0, [ak, 4'] - ^^'^ (k - k'). (17) 

There is some fiexibility in choosing the functions g{k,T) 
- or equivalently defining the states upon which the cre- 
ation and annihilation operators act. 

Specifically one can choose g{k,T) so that the Hamil- 
tonian operator commutes with the particle number op- 
erator, a]{j.ak, and all of the time dependence is carried by 
the field operator or one can choose g{k,T) so that the 
number of particles changes. The first option is tractable 
in closed form if Q is constant or piecewise constant in t; 
that is if TO 7^ and w = — 1 or m = and w is constant 
or piecewise constant. 



A. Vacuum Picture 

In the first case, the Hamiltonian operator is given by 

^y"d3fc fc(4ak + ak4) • (18) 

The Hamiltonian commutes with the particle number op- 
erator, so the universe remains in the vacuum state, |0). 
If the Hamiltonian takes this form, the functions g(fc,r) 
must satisfy the following differential equation 

%I> + (.»-2« ).(.,.,. (IS, 

If Q is constant we have 

g{k, t) = ^/7 [AJ,{kT) + BY,{kT)] (20) 



where v = ^/T+8Q/2. Two solutions of particular im- 
portance are for Q = 1 for w = 0, — 1 and ^ = to = and 
Q = for w = 1/3 and ^ = to = 0. We have 
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for = 1 Q. For w = — 1 (de Sitter space) one can use 
the general solution Eq. ((20|) when to 7^ and ^ 7^ 0, 
and during radiation domination {w — i), one can find 
a general solution in terms of Whittaker functions even 
when TO 7^ and f 7^ 0. In this paper the focus will be 
massless scalar fields, i.e. perturbations to the inflaton 
field and gravitational waves that can be modelled as 
two independent massless scalar fields [2]; therefore, the 
solutions given by Eq. ((2T|) and (|22|) will suffice for the 
radiation-dominated and de Sitter phases, respectively. 

In general if Q is only piecewise constant, one must 
insist that the function g(k, t) and its first derivative are 
continuous through the transition, because Eq. (fT9|l is 
second order. This determines the values of A and B in 
Eq. ((20)) after the transition {i.e. the Bogolubov trans- 
formation) . 

In particular if we look at the transition from vac- 
uum domination to radiation domination, we find that 
Eq. (j22p applies during vacuum domination and using 
the general solution, Eq (PO)) . we find that 
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during the subsequent radiation dominated epoch where 
a^max = l/(fc''RH) = 0-imH/k, the ratio of the Hubble 
parameter to physical wavenumber at the end of inflation. 

From Eq. (j23p we can read off the number of parti- 
cles that are created in a particular mode from an initial 
vacuum state; it is simply given by the squared of the co- 
efficient of the negative frequency term (see for example 
Eq. 1.26 of (H), 



Particle Picture 



(24) 



In this second technique one chooses a particular func- 
tion for g(k,T) and sticks with it for the entire calcula- 
tion. Specifically, we choose Eq. (pij) . This has two ad- 
vantages. First, this choice of g{k, r) is appropriate for 
radiation domination, i.e. the state of the universe after 
inflation - so the number of particles in a particular state 
remains constant after inflation. Second, this g{k, r) is 
appropriate for modes whose physical wavelengths are 
much smaller than the Hubble length, i.e. the situation 
in the distant past and the distant future (today). The 
particle number operator, aj^Ok, does not commute with 
the Hamiltonian during the de Sitter epoch, and its ex- 
pectation value gives the number of particles that we 
would expect to observe with a particle detector looking 
for excitations small compared to the present-day scale of 
the universe. Although there is considerable freedom in 
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choosing the functions g{k, t) due to the vacuum ambigu- 
ity that appears in quantum field theory in curved space- 
time, this particular choice of g{k,T) gives the canoni- 
cal quantization of a scalar field in a flat spacetime, in 
particular the locally flat, incrtial frame of a comoving 
detector, imparting some conceptual advantages of this 
picture over that standard picture outlined in §11 Al 

Because the Lagrangian for a scalar field in a 
Robertson- Walker spacetime evolves the same way as a 
scalar field in a flat spacetime but with a time-dependent 
mass, it is quite natural to express the Haniiltonian for 
the system as the sum of the flat space Hamiltonian [13] 
and time-dependent portion 



H = ^ 



Q 



:) (4< 



Q /. 



(25) 



Sk annihilates the vacuum of a radiation-dominated uni- 
verse, flk |0) = 0. 

The advantage of this technique is that it connects an 
initial state (the vacuum) to the final state in the basis 
in which we observe it. Because the final state is effec- 
tively the large-scale structure of the universe, following 
the quantum mechanical evolution of the scalar field be- 
comes exponentially cumbersome for modes much larger 
than the Hubble scale. Nevertheless, even the short-term 
evolution of modes as they expand beyond the Hubble 
scale provides new insights. 



C. Radiation-Dominated Spacetime Fock Space 

The Hamiltonian consists of a particle-conserving por- 
tion (which is the sum of the radiation-dominated result 
and a time-dependent correction) and a component which 
creates and destroys particles. The entire Hamiltonian 
commutes with the total momentum, so it is natural to 
examine how the Hamiltonian acts on states with zero 
total momentum. Let us take the following sum of states 



Ok' , Ok' 



^B„(t) |n, -k';n,k') , 



(26) 
(27) 



where we have assumed that the vacuum is normalised. 

|_B„(t)P is the probability of a comoving observer in 
the distant future detecting n particles with comoving 
momentum k' and n particles with momentum — k' if 
inflation ended at a conformal time r. This yields 



HI 
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r2A: 
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-2ifeT 
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+ (n + l)\n + 1, -k; n + 1, k) e^'''^ 



where Z = [dk, aj^] — (5(^'(k — k) is an infinite constant 
related to the vacuum energy of the system. Its value is 
independent of the number of particles in the states. 

We have replaced k' with k after performing the inte- 
gral over all momenta in Eq. (|25p . The integral with our 
choice of basis functions is manifestly even in k so the 
factor of one-half vanishes. Furthermore, the Hamilto- 
nian yields a solution separable in momentum-space, so 
we can solve for the evolution of each value of comoving 
momentum separately. As shown in the subsequent para- 
graphs, by solving this for a single value of the comoving 
momentum in de Sitter space one obtains solutions useful 
over a range of momenta. 

The time evolution of the coefficients Bnir) is given 
by the value of 



{n,-k;n,h\H\ib) 



.dB„(r) 
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+ B„+i(T)(n + 1)6-2*'=^ 

Since H evolves a state forward in time we obtain, 
.dAn{x) 



dx 



mA„_i(a;)e-2»7(:';)/:^ 



+ (n+l)A„+i(a;)e2^^(")/" 



(31) 



where we changed variables to a; = — l/(fcT) and absorbed 
an infinite phase into the definition of An, 

Ar,{x) = e-'(2"+z)(7(:'^)-i)/:^5^(a.) (32) 

7(t) = 2 + Qx^. (33) 

The variable x equals aH/k during the de Sitter phase 
and —aH/k during radiation domination. The value of x 
increases with time during both these phases. The phys- 
ical momentum p is given by p = fc/a, so p/H = l/\x\. 
The value of Q during the vacuum-energy-dominated 
phase {w = —1) is unity. During radiation domination 
{w — 1/3), in the massless case Q vanishes, so the coef- 
ficients An are constant during the radiation-dominated 
epoch (see Eg. [?T|) . 

It is important to keep in mind that the functional re- 
lationship between the effective mass and the conformal 
time dictated both the change of variables and the (in- 
finite) phase shift of the wavefunctions. In a situation 
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where Q is not constant with respect to t, a different 
set of substitutions are necessary. In general we must 
integrate up the time-dependence of the effective mass 
(Eq. I32p to rcnormahse the infinite zero-point energy of 
the field. This could be done numerically of course. 

Because ^„(r) differs from Bn only by a phase, 
|>l„(r)p still is the probability of a comoving observer 
detecting n particles with comoving momentum k and 
n particles with comoving momentum — k during radia- 
tion domination or the distant future if inflation ends at 
conformal time r. Furthermore, because the differential 
equation only depends on the product of fcr, the solu- 
tion to the equation gives the wavefunction in the Fock 
space of a radiation-dominated spacetime for a range of 
comoving momenta. 



D. Initial Conditions 

Initially, well before the mode crosses outside of the 
Hubble length, \Aq(t — oo)p = 1 and there are no 
particles. During inflation, the mode expands relative to 
the Hubble length. At the end of inflation, r = — trhj 
the universe reheats and the energy in the inflaton field 
decays into relativistic particles. After reheating we can 
no longer follow the quantum fluctuations of the inflaton 
field itself; however, any correlations in the inflaton field 
should also be present in the resulting perturbation to 
the density of relativistic matter. The quantum fluctua- 
tions of a fluid follow equations similar that of the scalar 
field [12], specifically the value of Q is the same for scalar 
fields and fiuids. We can also follow the evolution of the 
perturbations in a second scalar field which does not cou- 
ple to the inflaton as they cross back within the Hubble 
length. The perturbations induced by inflation on the 
metric itself can be decomposed into two minimally cou- 
pled scalar flelds , so the primordial spectrum of gravi- 
tational radiation falls into this category and furthermore 
these fluctuations remain in the linear regime until the 
present epoch. 

For simplicity, the vacuum and radiation-dominated 
regimes are pasted together so that after the end of in- 
flation r — trhj *-e- we assume that reheating is instan- 
taneous and that the universe does not expand during 
this time. The fact that inflation ends suddenly results 
in production of particles at all momenta. Specifically, 
the total energy of the particles produced diverges in the 
high-frequency limit; however, the sudden approximation 
is only valid for modes whose periods are much greater 
than the timescale for the transition. Higher frequency 
modes will evolve adiabatically so their final state will 
be the vacuum state. We also neglect any nonlinear cou- 
plings that may develop during reheating. This is proba- 
bly reasonable for tensor fiuctuations but non-linearities 
may play an important role for density fluctuations and 
this is an avenue for further work. 

The conformal time, r, continues to increase until the 
present day. The number of states occupied as well as 
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FIG. 1: Evolution of the expectation value of the particle 
number before and soon after Hubble length exit. The upper 
curve is for a four-particle initial state, the middle curve is 
a two-particle initial state and the lower curve is a vacuum 
initial state. 



the mean number of particles per mode increases dra- 
matically with the time that a particular mode spends 
outside of the Hubble length during an inflationary pe- 
riod. On the other hand, during radiation domination, 
the wavefunctions evolve only trivially, the coefficients 
An are constant. 



III. RESULTS 

As discussed in the previous paragraph, the natural 
initial condition for the simulation when a mode is much 
smaller than the Hubble length is = |0, k; 0, — k). We 
simulated the evolution of a massless real scalar field us- 
ing this initial condition and Eq. (|3ip . To assess the accu- 
racy of the simulations, the initial states 1 1, k; 1, — k) and 
|2, k; 2, — k) are also evolved forward. The inner product 
between these states and the initial vacuum state van- 
ishes initially. 



A. Number of Particles 

As Fig. [1] depicts, each initial state results in thousands 
of particles being produced over a wide range of comoving 
momenta. The rate of particle production is dramatically 
larger than that found by Mijic 9]. Also these calcula- 
tions do not exhibit the discontinuity that Mijic found at 
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H/p 

FIG. 2: Evolution of the scalar product between the various 
initial states. For clarity the value of 10 A^A„ — l) is 
depicted as well as the cross products. If the simulation were 
free of numerical errors, these values would all vanish. 

H/p = l/%/2 when the mass of a mode becomes imag- 
inary. Although the pair production is prodigious and 
the final states have little overlap with the initial states, 
the accumulated numerical error estimated by the scalar 
products and the conservation of the norm of the states 
is modest (Fig. [2]) . The accumulated error increases in 
a series of periodic steps rather than continuously. Each 
of these steps corresponds approximately to a conformal 
time where 2"/(x) = mr where n is a negative integer (see 
Eq.[33l). 

For all three sets of initial conditions, the number 
of particles increases dramatically for the modes that 
have exited through the Hubble length; however, in no 
case does the particle number increase exponentially with 
the comoving momentum; this would be indicated by a 
straight line in Fig. [TJ 

The simulations show that the occupation of the states 
is well described by a thermal distribution (Fig[3]), Such a 
distribution may be characterised by a single parameter, 
the expectation value of n. The total number of particles 
is of course N = 2n. Fig. |4] depicts the expectation value 
of 2n as a function of the ratio of the physical momentum 
to the Hubble parameter at the end of inflation. We have 
from Eq. (I24p . the number of particles in a particular 
mode 

Because the field theory for the field u is defined for a 
comoving volume, the physical energy density of particles 
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FIG. 3: The curves trace the final values of \A„\'^ as a func- 
tion of n/{n). The curves are normalised by the expectation 
value of n. Starting at the lower right corner of the panel and 
moving left, the curves trace modes with the following maxi- 
mum values of a; = H/p, 1.0, 1.5, 2.0 . . . 8.0. The distributions 
for X > 3 are nearly identical. 
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FIG. 4: The solid line traces the calculated value of (TV) = 
{2n} as a function the ratio of the physical momentum to 
the H is the Hubble parameter at the end of inflaton in the 
particle picture. The figure is similar to Fig. [T] but focuses on 
high frequency modes. The dashed line shows the value given 
by Eq. ([24)l in the vacuum picture. 
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produced within a range of physical momenta dp is 



dVdp 8tt^ p ^ ^ 

so the total energy density of the modes diverges as 
Inpmax/Pmin whcre Pmax ~ 1 / Ai is somc momentum cut- 
off, set by the duration of reheating [10|] . If the product 
of the physical momentum and the duration of reheating 
is much greater than one, the value of Q changes suf- 
ficiently slowly that the initial vacuum oscillation will 
evolve adiabatically to the final vacuum of the radia- 
tion dominated epoch. The infrared momentum cutoff 
is given approximately by the Hubble constant during 
reheating, p^^in ~ H ^lOi]. 

B. Primordial Correlations 

At some point the expansion of the universe becomes 
radiation dominated. If we assume that this transition 
is abrupt, we can calculate the subsequent evolution of 
a mode during the period of radiation domination as it 
shrinks relative to the Hubble length. The coefficients 
An remain constant during this time because Q = (see 
Eq. I3ip : in particular, the number of particles remains 
constant. When matter begins to dominate the energy 
density, the coefficients A„ again change with time as 
long as the mode lies outside the Hubble length because 
now Q ^ 0. The focus here will be modes that only spend 
at most a few Hubble times outside the Hubble length, 
so they will be well within the Hubble length at the onset 
of matter domination, and the values of An at reheating 
will characterise the fluctuations today. 

Even for modes that only spend a short time outside 
the Hubble length the final states long after reheating 
contain a large number of particles (Fig. [1]). However, 
the partition function for a given mode remains thermal. 
Fig. [3] shows that a simple rescaling of the distributions 
by the mean number of particles maps the final distribu- 
tions into each other. 

Although the partition function for the number of par- 
ticles is thermal for the final state of the fluctuations, the 
final wavefunctions for different values of the comoving 
momentum are correlated. 

The modes approach the classical regime when Xmax 
gets larger as Fig. [S] shows. The figure depicts the scalar 
product between the probability amplitudes of the final 
wavefunctions of modes that have spent nearly the same 
amount of time outside the Hubble length or have slightly 
different values of k. We find that for modes that only 
spend a short time outside the Hubble length the final 
state is quite similar over a wide range of momentum. 
The range over which the modes correlate decreases as 
2^max fo'' ^max ^1- As Xniax increascs the correlation 
between the modes diminishes. An important fact to 
keep in mind is that although the obvious correlations 
between the wavefunctions vanish as Xmax gets larger, 
the wavefunctions are still correlated albeit more subtly; 
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FIG. 5; The final value of for difi'erent values of 

— {H/p)niax- Starting with the uppermost curve on the 
left a;max = 1.6, 1.8, 2.0 .. . 4.4, 5, 6, 7, 8. 



the phases of the probability amplitudes are not strictly 
random but pseudo-random, that is to say calculable. 

Unfortunately, these correlated quantum phases are 
not directly observable. The wavefunction for each value 
of the comoving momentum resides in a separate Hilbert 
space, so our observations of the fluctuations on these 
small scales are drawn from the same underlying distri- 
bution but are otherwise independent. Only the most 
gross properties of the distribution are correlated such 
as the probability to flnd a certain number of particles 
in a particular state. On the other hand the correlated 
phases have a great pedagogical value in that they verify 
that the production of fluctuations during inflation does 
not require the destruction of information; it is strictly 
unitary. 



IV. SENSITIVITY TO THE DETAILS OF 
REHEATING 

Although the simulations here assumed that reheating 
quickly converts the potential energy of the scalar fleld to 
the kinetic energy of relativistic particles, simulations of 
a transition to a matter-dominated regime yielded very 
similar results. The flnal distributions are still thermal as 
in Fig.[3]and the correlations between the quantum states 
of modes with different values of k are similar to Fig. [5] 
except that in the case of a matter-dominated reheating, 
the correlations are typically weaker for given values of 
k and k'. 
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Of course, this is not a general demonstration that 
these results are robust with respect to the details of 
reheating but it does indicate that quantum correlations 
may persist to the present day. 

V. DISCUSSION 

The foregoing results show that fluctuations on scales 
of the present-day universe that passed through the Hub- 
ble length near the end of the inflationary epoch exhibit 
quantum mechanical correlations that may belie their 
birth from the vacuum. Today the comoving scale of 
the Hubble length at the end of inflation is 

i^'<To4v)"(T5S||v)"'™ <^«' 

where AI* is the vacuum energy associated with the in- 
flaton field and Trh is the temperature of the universe 
at the end of reheating. The analysis has assumed that 
reheating is quick and efficient [H, so Trh ~ M, 
yielding 

— = 2xl0^( ^ cm (37) 

The value of zrh for the comoving scale I/Zcrh is sim- 
ply one and k = fcRn/a^max- Consequently although the 
correlations are present on all scales, they are most ob- 
vious on the comoving scale of the Hubble length at the 
end of inflation {i.e. really small scales). On these small 
scales the density fluctuations are well into the non- linear 
regime today but tensor fluctuations, gravitational waves 
(GW), would still be a loyal tracer of these correlations. 

The expression given in Eq. [37] is very uncertain. A 
simple way to quantify the chance that we might observe 
these correlations directly is to calculate the number of 
e— foldings between the time when our present Hubble 
scale equalled the Hubble scale during inflation and the 
typical scale of future GW observatories probing inflation 
{e.g. the Big-Bang Observatory ^16;]): 

AN = ln^^^S8. (38) 



Typically today's Hubble scale is assumed to pass out 
through the Hubble length during inflation about 50—60 
e— foldings [l| before the end, so the millihertz scale 
would pass through the Hubble length 12—22 e— foldings 
before the end. However, the former number is highly 
uncertain. For example, if inflation occurs at a higher 
energy scale or if there is a epoch of late "thermal infla- 
ton" \rA [H, [l^, the number of e— foldings for today's 
Hubble scale could be as low as 25 

Even if we eventually observe the fluctuations on the 
comoving scales corresponding to the Hubble length at 
the end of reheating, one could argue that the wavefunc- 
tion only determines the probability of measuring a par- 
ticular amplitude and phase for the primordial gravita- 
tional waves at a particular momentum, and one would 
need to measure these fluctuations in several different 
realisations (universes) to unearth the correlations. 

The value of these correlations lies in that they verify 
that the production of fluctuations during inflation does 
not require the destruction of information; it is strictly 
unitary. Furthermore, one can use the wavefunction to 
calculate the density matrix for observations restricted to 
a portion of the spacetime {e.g. within the horizon |20|]) 
and estimate the entanglement entropy of the observable 
fluctuations. The calculations can also be easily extended 
to include non-linearities in the fleld either through grav- 
ity or self-interactions. Again these interactions generate 
entanglement entropy when fluctuations on a single scale 
are observed, giving the appearance of entropy produc- 
tion while the scale fleld remains in a pure state. Both 
of these calculations will be addressed in a future paper. 
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